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It is shown that a Bessel-like behaviour of the structure function F2 at small x, obtained 
for a flat initial condition in the DGLAP evolution equations, leads to good agreement 
with the deep inelastic scattering experimental data from HERA. 

1 Introduction 

The fairly reasonable agreement between HERA data PQ-[1] and the next-to- leading-order 
(NLO) approximation of perturbative QCD has been observed for Q 2 > 2 GeV 2 (see reviews in 
[5] and references therein) and, thus, perturbative QCD can describe the evolution of F 2 and 
its derivatives down to very low Q 2 values. 

The standard program to study the x behaviour of quarks and gluons is carried out compar- 
ing the experimental data with the numerical solution of the DGLAP equations [5] by fitting 
the QCD energy scale A and the parameters of the 2>profile of partons at some initial Qq 15] • 
However, to investigate exclusively the small- x region, there is the alternative of doing the 
simpler analysis by using some of the existing analytical solutions of DGLAP in the small- a; 
limit [S]-[T2] . It was pointed out in [5] that the HERA small- a; data can be well interpreted in 
terms of the so-called doubled asymptotic scaling (DAS) phenomenon related to the asymptotic 
behaviour of the DGLAP evolution discovered many years ago [12] ■ 

The study of |S] was extended in [TUJ-jT^] to include the finite parts of anomalous dimensions 
(ADs) of Wilson operators and Wilson coefficient^ This has led to predictions [TTJ fTJ] of 
the small- a; asymptotic form of parton distribution functions (PDFs) in the framework of the 
DGLAP dynamics, which were obtained starting at some Q 2 , with the flat function 

fc(Ql) = A a (hereafter a = q, g), (1) 

where f a are PDFs multiplied by x and A a are unknown parameters to be determined from the 
data. 

We refer to the approach of [10j-|12] as generalized DAS approximation. In this approach 
the flat initial conditions, Eq. ([lj, determine the basic role of the AD singular parts as in the 
standard DAS case, while the contribution from AD finite parts and from Wilson coefficients 
can be considered as corrections which are, however, important for better agreement with 
experimental data. 

* Speaker 

1 In the standard DAS approximation |13| only the AD singular parts were used. 



MPWLHC 2011 



1 



The use of the flat initial condition, given in Eq. ([I]), is supported by the actual experimental 
situation: low-Q 2 data [13 HE1 13] are well described for Q 2 < 0.4 GeV 2 by Regge theory with 
Pomeron intercept ap(0) = Xp + 1 = 1.08, closed to the adopted (ap(0) = 1) one. The small 
rise of HERA data [TJ 01 [THl HZ] at low Q 2 can be explained, for example, by contributions of 
higher twist operators (see [H]). 

The purpose of this paper is to demonstrate a good agreement [14] between the predictions 
of the generalized DAS approach [11] and the HERA experimental data [T] (see Fig. 1 below) 
for the structure function (SF) F 2 . We also compare the result of the slope 91nF 2 /91n(l/a;) 
calculation with the HI and ZEUS data [U |3J. Looking at the HI data [2] points shown in 
Fig. 2 one can conclude that X(Q 2 ) is independent on x within the experimental uncertainties 
for fixed Q 2 in the range x < 0.01. The rise of \{Q 2 ) linearly with InQ 2 could be treated 
in strong nonperturbative way (see [18] and references therein), i.e., X(Q 2 ) ~ l/a s (Q 2 ). The 
analysis [19j . however, demonstrated that this rise can be explained naturally in the framework 
of perturbative QCD. 

The ZEUS and HI Collaborations have also presented [3J the preliminary data for X(Q 2 ) 
at quite low values of Q 2 . The ZEUS value for X(Q 2 ) is consistent with a constant ~ 0.1 at 
Q 2 < 0.6 GeV 2 , as it is expected under the assumption of single soft Pomeron exchange within 
the framework of Regge phenomenology. It was important to extend the analysis of |19j to low 
Q 2 range with a help of well-known infrared modifications of the strong coupling constant. We 
used the "frozen" and analytic versions (see, [T4"]). 



2 Generalized DAS approach 



The flat initial condition ([!]) corresponds to the case when PDFs tend to some constant 
value at x — > and at some initial value Qq. The main ingredients of the results [TTJ[T2], are: 

• Both, the gluon and quark singlet densities are presented in terms of two components 
(" + " and " — " ) which are obtained from the analytic <5 2 -dependent expressions of the 
corresponding (" + " and " — ") PDF moments. 

• The twist-two part of the " — " component is constant at small x at any values of Q 2 , 
whereas the one of the " + " component grows at Q 2 > Qq as 



a = 2W d+s- D 



"/?o 



p 



In 



where a and p are the generalized Ball-Forte variables, 

P = as(Qo) ~ a s {Q ), 



la 'MQ 2 

Hereafter we use the notation a 

11 - 



12 



a 

21n(l/x) 
412 



27/3 
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(3) 



= a s /(47r). The first two coefficients of the QCD /3-function 
(2/3)/ and ft = 102 - (114/9)/ with / is being the number 



in the MS-scheme are 0o 
of active quark flavours. 

Note here that the perturbative coupling constant a s (Q 2 ) is different at the leading-order 
(LO) and NLO approximations. Hereafter we consider for simplicity only the LO approxima- 
tion^] where the variables a and p are given by Eq. Q when p = 0. 



2 The contribution of valence quarks is negligible at low x. 

3 The NLO results may be found in [TT] . 
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2.1 Parton distributions and the structure function F 2 

The SF F 2 and PDFs have the following form 

F 2 (x,Q 2 ) = e/,(x,g 2 ), f a (x,Q 2 ) = f+(x,Q 2 ) + f-(x,Q 2 ), (a = q,g) (4) 

where e = (%2{ e 2 )/ f is the average charge square. 
The small- x asymptotic results for PDFs f£ are 

f+(x,Q 2 ) = (A g + ^A q y (a)e-~ d ^ s + O(p), f+(x,Q 2 ) = LelM + ( p ), 

f;(x,Q 2 ) = - A - Aq e- d -^ s + 0(x), f~(x,Q 2 ) = A g e- d -^ s + 0(x), (5) 

where d-(l) = 16//(27/3 ) and d+(l) = 1 + 20//(27/3 ) is the regular part of AD d+{n) in the 
limit n — > lQ Here n is the variable in Mellin space. The functions l v (v — 0, 1) are related to 
the modified Bessel function l v and to the Bessel function J„ by: 

f ( \ — f if s > , v 

1,Aa) ~ \ i- v J v {ia), i 2 = -1, if s < • W 

2.2 Effective slopes 

As it has been shown in [TT], the behaviour of PDFs and F 2 given in the Bessel- like form by 
generalized DAS approach can mimic a power law shape over a limited region of x and Q 2 

/ a (z,Q 2 )~z~ Af(x ' Q2) and F 2 (z,Q 2 )~z~ A| * (X ' Q2) - 
The effective slopes X^(x,Q 2 ) and Ap ff (a;,(5 2 ) have the form: 



\ cff (r O 2 ) \ eS (r 2 \ - fsfaQ ) n H°) ~ n 1 



A9(X ' 0) " fq(x,Q 2 ) h(a)~ 41n(l/a;)' (?) 

where the symbol ~ marks the approximation obtained in the expansion of the modified Bessel 
functions, when the "— " component is negligible. These approximations are accurate only at 
very large a values (i.e. at very large Q 2 and/or very small a;). 



3 Comparison with experimental data 

Using the results of previous section we have analyzed HERA data for F 2 PQ and the slope 
d In F 2 / d \n(l / x) [HE] at small x from the HI and ZEUS Collaborations. In order to keep the 
analysis as simple as possible, we fix f — 4 and a s (M§) = 0.1166 (i.e., A^ 4 ) = 284 MeV) in 
agreement with the recent ZEUS results in pQ. 

4 We denote the singular and regular parts of a given quantity k(n) in the limit n — > 1 by k(n) and k(n), 
respectively. 
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Figure 1: x dependence of F2(x, Q 2 ) in bins of Q 2 . The experimental data from HI (open points) 
and ZEUS (solid points) [I] are compared with the NLO fits for Q 2 > 0.5 GeV 2 implemented 
with the canonical (solid lines), frozen (dot-dashed lines), and analytic (dashed lines) versions 
of the strong-coupling constant. 



As it is possible to see in Figs. 1 and 2, the twist-two approximation is reasonable at 
Q 2 > 2-^4 GeV 2 . At smaller Q 2 , some modification of the approximation should be considered. 
In Ref. |12j we have added the higher twist corrections. For renormalon model of higher twists, 
we have found a good agreement with experimental data at essentially lower Q 2 values: Q 2 > 0.5 
GeV 2 (see Figs. 4 and 5 in [12]). 

In Ref. [T3], to improve the agreement at small Q 2 values, we modified the QCD coupling 
constant. We consider two modifications. 

In one case, which is more phenomenological, we introduce freezing of the coupling constant 
by changing its argument Q 2 —> Q 2 + M 2 , where M p is the p- meson mass (see [14] and references 
therein). Thus, in the formulae of the Section 2 we should do the following replacement: 

a s {Q 2 ) -> a b (Q 2 ) ee a s (Q 2 + M 2 ) (8) 

The second possibility incorporates the Shirkov-Solovtsov idea [20J about analyticity of the 
coupling constant that leads to the additional its power dependence. Then, in the formulae of 
the previous section the coupling constant a s (Q 2 ) should be replaced as follows: (k = 1 and 2 
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Figure 2: As in Fig,l but for the Q 2 dependence of X e p 2 (x, Q 2 ) for an average small- a: value 
of x — 10~ 3 . The linear rise of X c p(x,Q 2 ) with InQ 2 [21 is indicated by the straight dashed 
line. For comparison, also the results obtained in the phenomenological models by Kaidalov et 
al. [35] (dash-dash-dotted line) and by Donnachie and Landshoff [53J (dot-dot-dashed line) are 
shown. 



at LO and NLO) 



where the symbol . . . stands for terms which are zero and negligible at Q > 1 GeV [50] at LO 
and NLO, respectively. 

Figure 2 shows the experimental data for X c p 2 (x,Q 2 ) at x ~ 10~ 3 , which represents an 
average of the x-values of HERA experimental data. The top dashed line represents the afore- 
mentioned linear rise of \{Q 2 ) with ln(Q 2 ). The Figs. 1 and 2 demonstrate that the theoretical 
description of the small-Q 2 ZEUS data for Xf 2 {x, Q 2 ) by NLO QCD is significantly improved by 
implementing the "frozen" and analytic coupling constants Uh(Q 2 ) and a an (Q 2 ), respectively, 
which in turn lead to very close results (see also |21j). 

Indeed, the fits for F 2 {x,Q 2 ) in O yielded Ql « 0.5-0.8 GeV 2 . So, initially we had 
X c p 2 {x,Qq) = 0, as suggested by Eq. (JIJ. The replacements of Eqs. (|8j) and ^ modify the 
value of X°p 2 (x,Qq). For the "frozen" and analytic coupling constants ct{ r (Q 2 ) and a an (Q 2 ), 
the value of Xp^(x, Qq) is non-zero and the slopes are quite close to the experimental data at 
Q 2 ss 0.5 GeV 22 . Nevertheless, for Q 2 < 0.5 GeV 2 , there is still some disagreement with the 
data, which needs additional investigation. 

For comparison, we display in Fig. 2 also the results obtained by Kaidalov et al. [52] and 
by Donnachie and Landshoff [23] adopting phenomenological models based on Regge theory. 
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Figure 3: The values of effective slope A§? as a function of Q 2 . The experimental points are 
same as on Fig. 4. The dashed line represents the fit from [2J. The solid curves represent the 
NLO fits with "frozen" coupling constant at x = 10 -2 and x = 10~ 5 . 



While they yield an improved description of the experimental data for Q 2 < 0.4 GeV 2 , the 
agreement generally worsens in the range 2 GeV 2 < Q 2 < 8 GeV 2 . 

The results of fits in [121 1 14] have an important property: they are very similar in LO 
and NLO approximations of perturbation theory. The similarity is related to the fact that the 
small-x asymptotics of the NLO corrections are usually large and negative (see, for example, 
a s -corrections [24] to BFKL approach [25] Then, the LO form ~ a s (Q 2 ) for some observable 
and the NLO one ~ a s (Q 2 )(l~ Ka s (Q 2 )) with a large value of K, are similar because A ^> Alc^] 
and, thus, a s (Q 2 ) at LO is considerably smaller then a s (Q 2 ) at NLO for HERA Q 2 values. 

In other words, performing some resummation procedure (such as Grunberg's effective- 
charge method 26J), one can see that the NLO form may be represented as ~ a s (Q 2 s ), where 
Q 2 S ^> Q 2 . Indeed, from different studies 27 , it is well known that at small-x values the 
effective argument of the coupling constant is higher then Q 2 . 

In the generalized DAS approach the small effect of the NLO corrections can be explained 
by separated contributions of the singular and regular AD parts. Indeed, the singular parts 
modify the argument of the Bessel functions (see Eq.Q) and the regular parts contribute to 
the front of Bessel functions [11] . 

Figure 3 shows the ^-dependence of the slope X e p(x, Q 2 ). One observes good agreement 
between the experimental data and the generalized DAS approach for a broad range of small-x 
values. The absence of a variation with x of Xp 2 (x,Q 2 ) at small Q 2 values is related to the 

5 It seems that it is a property of any processes in which gluons, but not quarks play a basic role. 
6 The equality of a s (Mg) at LO and NLO approximations, where Mz is the Z-boson mass, relates A and 
Alo: A (4) = 284 MeV (as in ZEUS paper on [l]) corresponds to A LO = 112 MeV (see [12]). 
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small values of the variable p there. 

From Figs. 2 and 6 in [TJ], one can see that HERA experimental data exists at x ~ 
1(T 4 -7- 10~ 5 for Q 2 = 4 GeV 2 and at x ~ 10~ 2 for Q 2 = 100 GeV 2 . Indeed, the correlations 
between x and Q 2 in the form x e g = a x 10~ 4 x Q 2 with a = 0.1 and 1 lead to a modification 
of the Q 2 evolution which starts to resemble InQ 2 , rather than lnlnQ 2 as is standard [T5] . 

4 Conclusions 

We have shown the Q 2 -dependence of the SF F 2 and the slope Xp = d\n F 2 /d\n(l/x) 
at small- a; values in the framework of perturbative QCD. Our twist-two results are in very 
good agreement with the precise HERA data at Q 2 > 2 GeV 2 , where the perturbative theory 
is applicable. The application of the "frozen" and analytic coupling constants af r (Q ) and 
«an(Q 2 ) improves the agreement for smaller Q 2 values, down to Q 2 > 0.5 GeV 2 . 

As a next step of investigations, we plan to fit the Hl&ZEUS data |3] and to extend the 
generalized DAS approach to evaluate the double PDFs which are very popular now (see [35] 
and references therein). Also we plan to use our approach to analyse the cross sections of 
processes studied at LHC by analogy with our investigations [2^] of the total cross section of 
ultrahigh- energy deep-inelastic neutrino-nucleon scattering. 

A.V.K. thanks the Organizing Committee of the 3rd International Workshop on Multiple 
Partonic Interactions at the LHC for invitation and support. This work was supported in part 
by RFBR grant 11-02-01454-a. 
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